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Several possible background sources determine the detectability of pep and CNO solar neutrinos 
in organic liquid scintillator detectors. Among such sources, the cosmogenic 11 C nuclide plays a 
central role. 11 C is produced underground in reactions induced by the residual cosmic muon flux. 
Experimental data available for the effective cross section for n C by muons indicate that 11 C will 
be the dominant source of background for the observation of pep and CNO neutrinos. 11 C decays 
are expected to total a rate 2.5 (20) times higher than the combined rate of pep and CNO neutrinos 
in Borexino (KamLAND) in the energy window preferred for the pep measurement, between 0.8 
and 1.3 MeV. 

This study examines the production mechanism of 1 C by muon-induced showers in organic liquid 
scintillators with a novel approach: for the first time, we perform a detailed ab initio calculation 
of the production of a cosmogenic nuclide, 11 C, taking into consideration all relevant production 
channels. Results of the calculation are compared with the effective cross sections measured by 
target experiments in muon beams. 

This paper also discusses a technique for reduction of background from 11 C in organic liquid 
scintillator detectors, which allows to identify on a one-by-one basis and remove from the data set 
a large fraction of 1 C decays. The background reduction technique hinges on an idea proposed by 
Martin Deutsch, who suggested that a neutron must be ejected in every interaction producing a 11 C 
nuclide from 12 C. 11 C events are tagged by a three-fold coincidence with the parent muon track and 
the subsequent neutron capture on protons. 
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I. INTRODUCTION 



Observation of solar neutrinos over the past 35 years 
by seven experiments has offered a unique opportunity 
to probe particle physics beyond the Standard Model of 
electroweak and strong interactions and physics of the 
stellar models. Radiochemical experiments have mea- 
sured the combined flux of a number of different neutrino 
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sources QH^EIOI- The only solar neutrinos targeted by a 
real time measurement thus far have been the 8 B neutri- 
nos above a detection threshold of about 5 MeV 0, E| 0] • 
Plans are in place to measure the sub-MeV neutrino spec- 
trum with organic liquid scintillator-based detectors, fo- 
cusing on the presumably abundant 7 Be neutrinos, with 
Borexino Igl, KamLAND and a possible detector at 
SNOLab [13. Such detectors also have the potential to 
probe the intermediate energy region, searching for the 
less abundant pep and CNO neutrinos. 

The possibility of detecting pep neutrinos in organic 
liquid scintillator based detectors is particularly intrigu- 
ing. It was recently pointed out that a measurement 
of the flux of pep solar neutrinos would yield essentially 
equivalent information about neutrino oscillation param- 
eters and the other solar neutrino fluxes as a measure- 
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ment of pp solar neutrinos at a comparable level of exper- 
imental uncertainty. Moreover, given the low theoretical 
uncertainty on the pep solar neutrinos flux, its measure- 
ment could allow investigating the matter-vacuum tran- 
sition region for solar neutrino oscillations . We recall 
that a matter-vacuum transition of solar neutrino oscil- 
lation is expected in the region between 2 and 3 MeV 
for the MSW-LMA solution of the solar neutrino prob- 
lem, and its observation would provide a further stringent 
test of the MSW-LMA solution. Due to this transition, 
the survival probability - i.e. the probability that elec- 
tron neutrinos emitted by the sun have not oscillated 
into other neutrino species when they arrive on earth - 
for pep neutrinos is about a factor of two larger than the 
one measured at higher energy by the SNO experiment 
for 8 B neutrinos 

The interaction rates for pep and CNO solar neutri- 
nos in organic scintillator predicted by the BP04 ver- 
sion of the Standard Solar Model by Bahcall and Pinson- 
neault [l2], in the currently preferred MSW-LMA solu- 
tion for the solar neutrino problem ^l]], are 2.1 and 6.6 
events per day in 100 metric tons of liquid scintillator 
respectively. When data recently released by the LUNA 
collaboration for the 14 N + p fusion cross section 0] are 
taken into consideration, the signal rate from CNO neu- 
trinos is significantly decreased to an expected 3.5 events 
per day in 100 tons. 

The pep neutrino energy spectrum is distinctive, with 
a single 1.44 MeV monochromatic line. The energy spec- 
trum for electrons scattered in v-e interactions presents a 
characteristic Compton-like edge at 1.22 MeV. Figure 
shows the expected (MSW-LMA) spectrum of scattered 
electrons from different neutrino sources in the energy 
range of interest for this study. We will focus our atten- 
tion on the detection and measurement of the pep neu- 
trino line. For this purpose we will set an observation 
window for the recoil electron between 0.8 and 1.3 MeV: 
for sake of simplicity, all event rates cited in the following 
will refer to this energy window, unless otherwise noted. 
The expected signal rate S in said window for pep and 
CNO neutrinos combined is 2.0 events per day in 100 
tons (1.5 events per day per 100 tons when using the 
most recent results from LUNA). 

There are three fundamental prerequisites for a succes- 
ful measurement of the pep and CNO solar neutrinos in 
an organic liquid scintillator detector. 

First, the internal background from long-lived radioac- 
tive sources must be carefully controlled. A 238 U and 
232 Th contamination at the 10~ 17 g/g level, coupled with 
a nat K contamination of 10~ 15 g/g would produce 0.6 
background events per day in 100 tons in the observa- 
tion window ; this is below the event rate expected 
for pep neutrinos. Contamination from long-lived radon 
daughters out of secular equilibrium with 238 U (in par- 
ticular from 210 Bi) must also be reduced below one count 
per day in the window of interest H, 0, 0] . 

A second prerequisite is a low external gamma-ray 
background from the construction materials of the de- 
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FIG. 1: Energy spectrum of electrons scattered by pep (dot- 
ted line) and CNO neutrinos (dash-dotted line) , for the MSW- 
LMA solution of the solar neutrino problem. The total spec- 
trum for all solar neutrinos (including pp, 7 Be, and 8 B neu- 
trinos, not shown separately) is also shown (continuous line). 
Neutrino fluxes for the spectra shown are from the BP04 ver- 
sion of the Standard Solar Model (see reference (Tljl l. The 
background signal expected from cosmogenic 11 C at Gran 
Sasso depth is superimposed (dashed line). For all the neu- 
trino and the C spectra, we assume that the energy reso- 
lution of the detector is 5% at 1 MeV and varies with the 
energy E as l/\/E. 



tector and from the surrounding rocks. The general 
strategy to solve the problem is to use pure buffer ma- 
terials to screen environmental radioactivity present in 
underground laboratories and extremely low radioac- 
tivity construction materials |l6l |. The typical spherical 
geometry of liquid scintillator based detectors has a pe- 
culiar effect on the spectrum of gamma-induced events 
reconstructed within a given radius: the spectrum gets 
harder at smaller radii. Therefore, in the innermost 
part of the detector the 7 ray-induced background in 
the pep window [0.8-1.3 MeV] is larger than the corre- 
sponding background in the 7 Be window [0.25-0.8 MeV]. 
This potential problem may be counterbalanced by re- 
defining the fiducial mass for the observation of pep and 
CNO neutrinos. For reference, in the 100 tons fiducial 
mass of Borexino for Be neutrinos detection the exter- 
nal background in the [0.8-1.3 MeV] window is expected 
at 1 event/day compared with a neutrino signal of 

2 events / day. Reducing the mass to 70 tons would lower 
the background by a factor 10 0, ^| while losing only 
30% of the signal. 

The third fundamental condition for the observation 
of pep and CNO neutrinos is a low internal cosmogenic 
background production. This topic is the main subject of 
this paper. Section [H] is an introduction to the problem 
of 11 C cosmogenic background in liquid scintillator de- 
tectors. Section lnTl describes the production channels for 
cosmogenic 11 C and the results of our calculations for the 
expected production rate in a 100 tons liquid scintillator 
target. Section llVl offers a comparison of our estimate 
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with another estimates available in the literature for the 
production of 11 C by direct interaction of a muon with a 
12 C nucleus through virtual photons. Finally, section Ivl 
proposes a veto mechanism for 11 C production events, 
based on the double coincidence with the parent muon 
and a neutron produced in the reaction. We draw our 
conclusions in section IvTl 



II. COSMOGENIC BACKGROUND IN DEEP 
UNDERGROUND DETECTORS 

Cosmogenic radioactive nuclides are produced in deep 
underground detectors in reactions triggered by the resid- 
ual muon flux. As originally pointed out in [T3 , the fun- 
damental scale of the process is given by the neutron pro- 
duction rate: neutrons are an important and ubiquitous 
by-product of cosmic ray-induced nuclear reactions. 

A comprehensive review of the experimental results 
on neutron production in underground laboratories can 
be found in |J^. The first measurement at Gran Sasso 
(depth of 3800 meters of water equivalent, muon flux 
of 1.2 m -2 hr _1 , average muon energy 320 GeV 0|) 
was performed with the CTF detector [l5j in 1995 yield- 
ing (1.5 ± 0.1) x 10~ 2 neutrons//x/m (neutrons per me- 
ter of fj, track) in scintillator with a density of 0.88 
g/cm 3 H El- A later measurement by the LVD 
experiment gave (1.4 ± 0.4) x 10~ 2 neutrons/^t/m [22^ . 
The figure from the CTF experiment translates into a 
neutron capture rate of 40 events per day in 100 tons of 
scintillator. 

Deutsch suggested |2(| that the underground produc- 
tion rate for all of the most significant cosmogenic nu- 
clides out of a target mass composed of 12 C and 1 H could 
be estimated from the neutron production rate alone. 
The list includes 8 Li, 9 Li, n Bc, 8 B, 12 B, and 9 C, all 
with mean lives below one minute and easily taggable 
with the parent muon. 7 Be is the cosmogenic radionu- 
clide with the longest mean life, 77 days. n C also poses 
problems, given its 30 minutes mean life which does not 
allow identifying its decays by tagging them with the par- 
ent muon alone. Deutsch hypothized that 11 C would be 
one of the most likely cosmogenic nuclides by-products 
of muon-induced cascades, and estimated that 11 C would 
account for 5% of the total neutron production rate. This 
would correspond to 2 events per day in the 100 fiducial 
tons of Borexino. 

The inclusive cross section for the production of sev- 
eral cosmogenic nuclides in muon-induced cascades was 
measured by a target experiment on a muon beam at 
CERN . The experiment used a liquid scintillator 
target. The muon showers were built up in 240 cm of con- 
crete and 200 cm of water used as absorbers and placed 
in front of the 12 C targets. The use of two positively 
charged muon beams of 100 and 190 GeV allowed to 
extract information about the energy dependence of the 
inclusive cross sections. The results in |23j have shown 
that 7 Be is not among the most likely products of cos- 



TABLE I: Depth, residual muon flux, average muon energy, 
and neutrons capture rate (N) at KamLAND, Borexino, and 
SNOLab. Data on muons are from references p. ll9ll2Ml24ll25l 
I26l . l27j| . The neutron capture rate for Borexino comes from the 
CTF experiment HIHEl 

and the one for KamLAND is the 
rate measured in the detector and reported in reference 
The neutron capture rate for SNO is by extrapolating the 
CTF data with the scaling law introduced in 1171 ; when the 
same procedure is applied to calculate the capture rate in 
KamLAND, the result of 294 cts/d/100 tons is fully consistent 
with the value measured by KamLAND. 
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2700 


9.6 


285 


300 


Borexino 


3800 


1.2 


320 


40 


@ SNOLab 


6000 


0.012 


350 


0.43 



TABLE II: Total expected 11 C decay rate and "C-induced 
raw background rate (Bo) in the pep window [0.8-1.3 MeV] 
at KamLAND, Borexino, and SNOLab (approximately 35% 
of 11 C decays produce an event the pep window). Data for 
Borexino and KamLAND are from reference |23j|. extrapo- 
lated from data of target experiment on muon beam at 100 
and 190 GeV. Data for SNOLab are calculated from the muon 
flux and average energy reported in table[TJusing the same ex- 
trapolation method of reference [23l | . 
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KamLAND 


107 48 


37 


Borexino 


15 52 


5.1 


@ SNOLab 


0.15 55 


0.056 



mic ray-induced reactions: its production rate at Gran 
Sasso depth is expected to be less than 0.1 events per 
day in 100 tons I23. A precise measurement of the in- 
clusive n C production rate was also made available. At 
the mean energy of muons at Gran Sasso, the production 
rate is 15 events per day in 100 tons |2^] or, equivalently, 
52 x 10~ 4 "C/^/m. 

The rate is different for other locations due to the dif- 
ferent muon rates and muon average energies. Muon 
rates and expected 11 C production rates in KamLAND, 
Borexino, and at SNOLab are summarized in tables [lj 
andUD 

III. 11 C PRODUCTION IN MUON-INDUCED 
SHOWERS 

n C is a positron emitter with a 0.96 MeV end point; 
approximately 35% of its decays produce an event in the 
pep window [0.8-1.3 MeV]. The raw background rate (B ) 
from 11 C in the pep window varies with the depth of the 
location of the experiments. Values for the raw back- 
ground for Borexino, KamLAND, and at SNOLab are 
also reported in table ITU The background at SNOLab is 
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FIG. 2: Cross sections for n C production from 12 C as a 
function of energy. 

sufficiently low to enable pep and CNO neutrinos obser- 
vation without need of any cuts on 11 C. 

Deutsch pointed out that the only way to create 11 C 
is to knock a neutron off the 12 C nucleus and suggested 
to look for a neutron in the final stated of the reaction, 
emphasizing the possibility of a three-fold coincidence 
with the parent muon track and the neutron capture on 
protons in the scintillator to tag the 11 C events on a 
one-by-one basis [2(|. On the other hand, as suggested 
by Calaprice, there is also the possibility of creating a 
n C while ejecting a deuteron in a (p,d) exchange reac- 
tion. The (p,d) reaction would create an invisible channel 
for 11 C production since the nuclide produced through 
such process cannot be tagged with the three-fold coinci- 
dence mentioned above. Similarly, reactions triggered by 
7T mesons can also produce decays in invisible channels, 
as explained later. 

Following is a list of the leading reactions that can pro- 
duce 11 C, together with their cross section and references 
to specific studies. For all the reactions with a neutron 
in the final state, the energy threshold is ~ 20 MeV (i.e. 
the neutron binding energy in 12 C). 

12 C(7,n) 11 C |2j|: the cross section for the 7-ray in- 
duced process peaks at 7 mb around 23 MeV. The 
cross section value in the peak region, relatively 
large for an electromagnetic interaction, is due to 
the nuclear giant dipole resonance |2£j . 

12 C(n,2n) 11 C 30]: the cross section has a sharp peak 
(17 mb) around 33 MeV. At higher energies, we rely 
on the set of experimental data by Kim et al. |30j in 
the range 40-150 MeV, which is the only available 
for energies above 40 MeV. Those data are affected 
by a large experimental error, as high as 40%. As 
pointed out by the authors, their results disagree 
starkly with theoretical expectations which predict 
much lower values for the cross section in the same 
energy range Also, we assume that at energies 
higher than 150 MeV the cross section keeps the 



constant value attained in the range 70-150 MeV. 
Uncertainties in the knowledge of the cross section 
for this process represent the largest systematic er- 
ror in our ah initio calculation of the rate of pro- 
duction of cosmogenic 11 C. We estimate that the 
systematic error attributable to this source could 
reach 50% of the production rate expected from 
the 12 C(n,2n) 11 C channel. All the other cross sec- 
tions are known with a precision better than a few 
percent. 

C(p,p+n) n C |33|: the cross section reaches a peak 
value of 98 mb at 40 MeV, and decreases to a 
plateau of 30 mb, constant up 1 GeV. 

C(p,d) 11 C [33j]: the cross section has a threshold of 
16 MeV and has been measured at 52 and 65 MeV 
(15 and 10 mb respectively). The only measure- 
ment available in the range above 100 MeV tells us 
that the cross section is in the range of a few /xb, 
and therefore negligible. 

C(7r _ ,7r _ +n) 11 C [34]: the cross sections exhibits a 
broad peak centered around the (3,3) resonance 
for the pion-nucleon quasi-elastic interaction (see 
Dropesky et al. [34|]) with a value of 70 mb at 
190 MeV. Data are available up to 550 MeV, and 
show that the cross section reaches a plateau above 
400 MeV. We assume that the cross section keeps 
a constant value at higher energies. 

C(7r + ,7r+N) 11 C [34(: the cross section exhibits a 
broad resonance peak in the same region, reaching 
45 mb around 160 MeV. Data are available up to 
470 MeV, and show that the cross section reaches a 
plateau above 350 MeV. We assume that the cross 
section keeps a constant value at higher energies. 
Contrarily to n~ , positive mesons do not necessar- 
ily produce a neutron in the final state: the N in 
the final case stands for nucleon and can be either 
a proton or a neutron. Due to the possibile charge 
exheange occurring in the strong meson-nucleon in- 
teraction, the fragments in the final state can be ei- 
ther (7r + -|-n) or (7r°+p), the latter having a thresh- 
old of 13 MeV. No data are available on the relative 
composition of the final state, but it is expected 
from theoretical considerations (see Chivers et al. 
in [34| ) that the frequency of the invisible channel, 
with a proton in the final state, should account for 
2/3 of all the 7r + -induced reactions. We rely on this 
assumption in our calculations. 

C(e,e+n) 11 C [3^: The direct interaction of electrons 
and positrons with a nucleus, through a virtual 
photon, is expected to have a small cross section, 
of the order of a (= 1/137) times the peak value of 
the cross section for real photons. The measured 
value of the cross section at 30 MeV is 15 /ib. At 
higher energies, the cross section can be calculated 
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by using the von Weizsacker approximation 36] : 

o c = J N(v)a^(v)du (1) 

where a c is the cross section for the 11 C produc- 
tion induced by electrons, er 7 is the cross section 
for n C production by real photons, v is the energy 
of virtual photons and N is the number of virtual 
photons. In case of high energy, ultra-relativistic 
charged particles inducing nuclear reactions with 
low momentum transfer, the number of virtual pho- 
tons can be approximated by [37| : 

N(z/) = (a/irv) [2 In (E/m) - 1] (2) 

where E is the energy of the charged particle and 
m is its mass. The cross section for the production 
induced by electrons is then given by: 

a c = (a/n) [2 In (E/m) - 1] <r_i (3) 

where <r_i = J dvo^iy) jv is the inverse energy 
weighted moment of the photodisintegration cross 
section er 7 . Using the photodisintegration cross sec- 
tions from references |28|, we calculated a value of 
1.74 mb for With this value, the cross sections 
for the production of n C by electrons of 100 MeV, 
1 GeV, and 10 GeV are 39 fib, 57 ,ub, and 76 ^b 
respectively. The values of the cross sections at in- 
termediate energies are interpolated from the above 
values. 

12 C(/x,/i,+n) 11 C The direct interaction of muons with 
a nucleus, through a virtual photon, is usually re- 
ferred to as "muon spallation" . The cross section 
for the process can be calculated using the same 
procedure detailed above for electrons. The result 
is a cross section of 58 ^b for muons at 320 GeV. 
In the range 100-350 GeV the cross sections have 
values very close to the one just quoted, given the 
logarithmic dependence of the number of virtual 
photons N upon the energy E, as shown in equa- 
tion 121 

The cross sections for 11 C production by photons and 
hadrons, compiled using the procedure detailed above, 
are shown in figure |21 

We performed a full simulation of muon-induccd show- 
ers with the particle transport code FLUKA 38]. The 
FLUKA code has been used by Wang et al. [l8] to cal- 
culate the production rate of neutrons by muons in liq- 
uid scintillator at several depths, and has been found 
to reproduce experimental results very well. Recently, 
FLUKA has been used by Kudryavtsev et al. [27| to cal- 
culate the distance between the parent muon track and 
the point of capture on protons of neutrons produced 
in scintillator by muon-induced cascades at Gran Sasso 
depth, and results were found to be in agreement with 
the experimental data from the LVD experiment |22|. 
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FIG. 3: Cumulative range of secondaries generated in showers 
induced by negatively charged muons at 320 GeV. Results are 
quoted in cm of range per meter of fi track. 



We used FLUKA to calculate production rates and 
ranges of all the prominent secondaries, i.e. protons, neu- 
trons, 7r mesons, and 7-rays. We simulated showers orig- 
inating from negatively charges muons (fJ.~) at 100 and 
190 GeV (the energies of the muon beams for the ex- 
periement described in (23j), at 285 GeV (average energy 
at Kamioka), at 320 GeV (average energy at Gran Sasso), 
and at 350 GeV (average energy at SNOLab) . The target 
material in the simulation was the solvent of the liquid 
scintillator for Borexino, trimethylbenzene (C9H12), with 
density 0.88 g/cm 3 (incidentally, this makes up 20% of 
the solvent used in KamLAND Q). Results should not 
vary greatly with other organic solvents, given that typi- 
cal values of density and mass ratio between carbon and 
hydrogen are close to the values of trimethylbenzene. We 
tracked muons for 100 meters, and for each of the promi- 
nent secondaries we calculated the cumulative range of 
the particles as a function of the particle energy with a 
10 GeV cutoff. The results for secondary particles with 
energy below 1 GeV in showers induced by negatively 
charged muons at 320 GeV are shown in figure [3] 

We then turned to the computation of the 11 C produc- 
tion rate from each one of the interactions listed above, 
for each of the energies of the muons taken into con- 
sideration. The production rate has been calculated by 
taking the energy convolution, for each of the possible 
interactions, of the cross sections with the range of the 
secondary particles responsible for inducing that partic- 
ular interaction. Results are summarized in table II I II 
The production rate in the invisible channels rate has 
been calculated by adding the rate from the 12 C(p,d) 11 C 
reaction to 2/3 of the rate from the 12 C(7r+,7r+N) 11 C 
reaction. 

The error quoted in table Hill accounts for the system- 
atic error, coming from two main sources. The most im- 
portant source is the uncertainty in the knowledge of the 
cross section for the process 12 C(n,2n) 11 C. Our best es- 
timate for the systematic error is 50% of the production 
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TABLE III: Production rates for 11 C in muon induced show- 
ers. The calculated total production rates are compared with 
the experimental values available at 100 and 190 GeV from 
Ref. |23| . and with the extrapolated values at the mean muon 
energy for KamLAND and Borexino (also from Ref. |23|) and 
at SNOLab. The procedures used to determine the expected 
rate from the invisible channels and the systematic error af- 
fecting our calculation are outlined in the text. 
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12 C(p,d) 11 C 


0.2 


0.4 


0.5 


0.6 


0.6 


12 C( 7 ,n) 11 C 


19.3 


26.3 


33.3 


35.6 


37.4 


12 C(n,2n) n C 


2.6 


4.7 


7.0 


8.0 


8.2 


12 C(7r+ ! 7r+N) 11 C 


1.0 


1.8 


2.8 


3.2 


3.3 


12 C(7r-,7r-+n) 11 C 


1.3 


2.3 


3.6 


4.1 


4.2 


12 C(e,e+n) n C 


0.2 


0.3 


0.4 


0.4 


0.4 


12 C(^+n) n C 


2.0 


2.3 


2.4 


2.4 


2.4 


Invisible channels 


0.9 


1.6 


2.4 


2.7 


2.8 


Total 


28.3 


41.3 


54.8 


59.9 


62.2 


la systematic 


1.9 


3.1 


4.4 


5.0 


5.2 


Measured 

lcr experimental 

Extrapolated 


22.9 
1.8 


36.0 
2.3 


47.8 


51.8 


55.1 



rate in such channel, which is between 5% and 7% of the 
total production rate depending on the energy. We also 
make the conservative assumption that the systematic 
error associated with the cross sections used by FLUKA 
to calculate the range of the secondaries accounts for 5% 
of the total production rate. We combine the two as 
independent sources of error. The statistical error asso- 
ciated with our Monte Carlo calculation is 0.6% and is 
negligible with respect to the systematic error. 

The total calculated rate is systematically 20% higher 
than the measured rates on a beam at 100 and 190 GeV, 
even though the values are still within twice the combined 
experimental and systematic uncertainties. A possible 
explanation for the systematic discrepancy is the differ- 
ence in the geometry between the experiment and the 
simulation (in the experiment, the muon shower is built 
up in 240 cm of concrete and 200 cm of water placed in 
front of the scintillator target; in our simulation we con- 
sider a bulk scintillator volume). The calculated rates at 
285, 320, and 350 GeV, although systematically higher 
than the extrapolated values, are still in good agreement 
with them when considering the systematic error. 

We also calculated the production rates for positively 
charged muons (/i + ) at 320 GeV. The difference in the 
production of 11 C, for all of the channels considered, by 
/i + and fj," at 320 GeV is within the statistical error of 
our Monte Carlo calculation. We conclude that the de- 
pendence of the inclusive cross section for 11 C production 
on the sign of the charge of the muons is negligible within 
the scope of the study presented in this paper. 

We took into account also the possibility of producing 



n C nuclides from the target nuclide 13 C, which has a nat- 
ural isotopic abundance of 1.1% 39). The contributions 
of channels such as 13 C(7,2n) n C 40j, 13 C(7r+,d) 11 C [Ilj, 
and 13 C(p,t) 11 C (4^, are negligible with respect to the 
corresponding rates for processes with the same incident 
particles on the target nuclide 12 C, owing to smaller cross 
sections and to the low natural isotopic abundance of 13 C. 

As shown in figure |3 the particle content of muon- 
induced showers is, as expected 01, dominated by 7-rays 
and electrons. Around 25 MeV, at the giant dipole res- 
onance of 12 C where the cross section for 7-rays is quite 
large, the cumulative range of 7-rays is two to three or- 
ders of magnitude larger than the corresponding values 
for neutrons, and four to five orders of magnitude larger 
than for other hadrons. As a consequence, the domi- 
nant process for the production of 11 C nuclides is the 
(7,n) exchange reaction, accounting for ^60% of the to- 
tal production rate, even if the hadronic cross sections for 
the 11 C production are up to a factor 10 larger than the 
peak value for 7-rays. Electrons are not as effective as 
7-rays because their cross section is lower by a factor a. 
Other hadronic channels with a neutron in the final state 
and the muon-induced photodisintegration, also carrying 
a neutron, account for an additional 35% of the total 
production rate. The rate of 11 C production in the two 
invisible channels, corresponding to the (p,d) and the 
(7r + ,7r°-|-p) exhange reactions, accounts for only about 
5% of the total production rate. 

Deutsch's idea of eliminating the 11 C events by look- 
ing at the three-fold coincidence turns out to be still 
valid even in presence of invisible channels: only 1 out 
of 20 11 C nuclides is produced without a neutron in the 
final state. In section[V]we will quantify the effectiveness 
of the 11 C background reduction technique based on the 
double muon+neutron tagging. 



IV. COMPARISON WITH PREVIOUS 
ESTIMATES 

O'Connell and Schima |44| calculated the production 
rate of several radioactive nuclides in carbon, oxygen, 
and argon targets, at sea level and at the depth of Kam- 
LAND, only for the photoproduction induced by the vir- 
tual photons associated with the muons. For the n C 
nuclides, they took into consideration only the chan- 
nel 12 C(/i,/i+n) n C which results, in our estimate, in 
a 5% contribution to the total production rate. They 
estimated a production rate of 11 C of 15 events per 
day in 100 tons of carbon at the KamLAND depth, 
or, scaling with the muon flux quoted in table [IJ about 
6.0 x 10~ 4 nuclides//z/m track in organic liquid scintil- 
lators. Their result is to be compared with our estimate 
for the same channel reported in tablc lTTTl which amounts 
to 2.4 x 10~ 4 nuclides/^t/m. The estimate of O'Connell 
and Schima is about 2.5 times higher than our estimate. 
The reason for the discrepancy lies in the value chosen 
for (7-1, the inverse energy weighted moment of the pho- 
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todisintegration cross section. The value of a_i quoted 
in m[ is 4.5 mb, while the value we calculated from the 
cross sections in use in this paper is 1.74 mb. It is worth 
noting that work cited as the source for the photoneutron 
cross sections in the paper of O'Connell and Schima, the 
Atlas of photoneutron cross sections obtained with mo- 
noenergetic photons |45| . offers two values for o~-\ extrap- 
olated from three different experiments, whose average is 
1.65 mb. This is in excellent agreement with our calcula- 
tion and is about a factor 2.5 lower than the value quoted 
in 0. 



V. 11 C REDUCTION TECHNIQUES 

As shown in the previous section, at least one neutron 
is produced in association with 95% of the n C nuclides. 
This fraction of the 11 C background can be lowered with 
the muon+neutron tagging. The remaining 5% cannot 
be reduced with this technique. 

In order to suppress the 11 C background, one needs to 
identify the position and time of each neutron created by 
a muon-induced shower and then captured on protons in 
the scintillator. Neutron capture on protons results in the 
emission of a distinctive 2.2 MeV 7-ray. In hydrocarbons, 
neutrons can also be captured on 12 C resulting in 7-rays 
of combined energy 4.9 MeV |4fil]. The cross section for 
the capture on 12 C is ~1% 47] of the cross section for 
capture on protons. Once the neutron capture time and 
position are known, one needs to apply a cut in space 
and time around every capture point. The events within 
a time t from the double muon+neutron coincidence and 
inside a sphere of radius r from the neutron capture point 
are rejected. This technique, originally suggested in [2(| , 
has recently been succesfully applied in the 4-ton proto- 
type Counting Test Facility of Borexino . 

The length of time t for which events are rejected 
should be set to a few times the mean life of 11 C. Note 
that the information carried by the neutron capture does 
not tell us about the position of the 11 C birthplace. 
Therefore, it is important to set the radius r of the spher- 
ical cut to a few times the average neutron range (Note: 
the spatial resolution of organic liquid scintillator detec- 
tors is about 10 cm for the 2.2 MeV 7-rays from neu- 
tron capture and can be neglected with respect to the 
average neutron range). For this reason, we calculated 
with Monte Carlo methods the energy distribution of the 
neutrons produced in association with 11 C nuclides for 
Borexino, KamLAND and at SNOLab. The calculation 
was performed using the relative weights of the 11 C pro- 
duction in the different channels as determined for muons 
of 320 GeV. The procedure is correct for all the three lo- 
cations since the relative weights 11 C production in the 
different channels at 285 and 350 GeV are within 3% of 
the values at 320 GeV (see table llllf) . We used the re- 
sulting distribution, shown in figure^ to source neutrons 
into FLUKA, calculating the distribution of the range of 
neutrons produced in association with 11 C, shown in fig- 
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FIG. 4: Energy distribution for neutrons produced in asso- 
ciation with 11 C nuclides. 
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FIG. 5: Range of neutrons produced in association with C 
nuclides. 

ure[3 

The average energy of neutrons produced in associa- 
tion with 11 C is much lower than the average energy of all 
neutrons produced in muon-induced showers: this is due 
to the dominant production mechanism being the pho- 
toneutron reaction at the giant dipole resonance of 12 C 
at 23 MeV. Correspondingly, the average range of neu- 
trons associated with 11 C production is also much lower 
than the average range of all the neutrons produced in 
the shower. 

In our calculation we assumed that the 11 C nuclides 
displacement between the point where it is created and 
the point where it decays is negligible with respect to the 
range of neutrons created in association with 11 C nuclide. 
For that to happen, the convective motion of the scintil- 
lator has to be sufficiently slow. This can be achieved, 
for example, by maintaining a small temperature gra- 
dient pointing upward everywhere in the detector. The 
average range of neutrons created in association with n C 
nuclides, whose distribution is shown in hgure[Sl is 44 cm. 
KamLAND data 49] show that the measured average dis- 
placement of the diffusive 222 Rn over its 5.5 days mean 
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life is less than 1 meter. Therefore the assumption that 
the 11 C nuclides displacement over their 30 minutes mean 
life can be kept small with respect to the neutrons range 
seems fully justified. 

Due to the presence of the (n,2n) exchange reaction 
that yields two neutrons in the final state, an average 
of 1.14 neutrons are created in interactions producing 
n C nuclides. Those neutrons, if sufficiently energetic, 
can also trigger nuclear reactions knocking off other neu- 
trons: our calculation indicates that an average number 
of about 1.2 neutrons are captured in the scintillator for 
each neutron produced in a C-forrrring reaction. 

The efficiency e of rejecting 11 C events tagged with 
the muon+neutron coincidence is equal to the combined 
efficiencies for the cut in space, C( r )i an d in time, rj(t): 

e = C(r) ■ ry(t) (4) 

where r](t) is the efficiency for the rejection of 11 C events 
when a cut in time for a time span equal to t is applied 
around the neutron capture point. Given the expected 
decay distribution of n C, we obtain: 
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FIG. 6: Correlation between the signal to reduced background 
ratio (R) and the dead mass-time fraction for Borexino, as a 
function of spatial (£) and time (??) efficiencies of 11 C events 
rejection. The only background considered in the pep energy 
window [0.8-1.3 MeV] is the 11 C background surviving the 
cuts described in the text. Isocontours for the ratio R (for 
values R = 1, 2, 5) and for the dead mass-time fraction D (for 
values D = 1%, 10%, 50%) are shown in the graph. 



-t/T 



(5) 



where r is the mean life of 11 C, 30 minutes. 

In equation £(r) is the rejection efficiency for n C 
when a cut in space is applied around the neutron cap- 
ture point, corresponding to a sphere of radius r centered 
around the capture point position. Given the distribu- 
tion n(s) for the range of neutrons shown in figure 03 the 
efficiency for the rejection is equal to: 



C(r) 



n(s)ds 



n(s)ds 



(6) 



In order to quantify the effectiveness of the background 
reduction technique we introduce a figure of merit R, in- 
dependent from the experiments and their locations. R is 
defined as the ratio of the pep+CNO neutrino signal rate 
(S) to the residual background rate (B) from n C after 
suppression of all 11 C events identified through the the 
muon+neutron coincidence (both S and B are computed 
in the pep energy window [0.8-1.3 MeV]). Note that the 
figure of merit R accounts only for the remaining back- 
ground from 11 C, but other possible and independent 
sources of background are neglected in its definition. The 
general expression for R is: 



R = S/B 



S/B 



F + (1 - F) (1 — (rj) 



(7) 



where F ~ 0.05 is the fraction of n C production rate 
in invisible channels; S (signal rate) was discussed in sec- 
tion^and B (raw n C background rate, experiment- and 
location-dependent) was discussed in section IIIII and in 



table [n] Formula? 0] and [7\ show that there is a one-to- 
one correspondance between the figure of merit R and 
the combined rejection efficiency e. 

Using equation the signal rate reported in section [I] 
and background rates from table[n] we estimate that the 
optimal signal to background ratio achievable in Borex- 
ino (KamLAND) is 8 (1) in case all of the n C associated 
with a neutron are successfully tagged. The optimal limit 
cannot however be reached due to dead mass-time limi- 
tations, as shown below. 

We define the dead mass-time fraction D as the fraction 
of (mass x time) data taking lost to the space and time 
cuts around a neutron capture event. Treating rj and 
£ as independent variables, t and r can be derived by 
inverting equations [S] and El and become functions of rj 
and £ respectively. We can then compute the dead mass- 
time fraction D corresponding to the chosen values of the 
cuts £ and rj as: 



D = 1 



(8) 



where p is the scintillator density, and N is the neutron 
rate per unit mass in the detector. The expected neutron 
rates for the experiments taken into consideration are 
summarized in tabled 

Figure shows the correlation between the signal to 
reduced background ratio R and the dead mass-time frac- 
tion D for Borexino, as a function of spatial (£) and 
time (77) efficiencies for n C events rejection. In table HVl 
we show the optimal values of the dead mass-time frac- 
tion D for the three experiments for fixed and given val- 
ues of the ratio R between signal and reduced 11 C back- 
ground. The values are calculated by minimizing the 
dead mass-time fraction D while keeping the product of 
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TABLE IV: Optimal values of the dead mass-time fraction 
(D) for KamLAND, Borexino and SNO as a function of the 
figure of merit R, defined as the ratio between the signal rate 
(S) of pepand CNO neutrinos (BP04 model for MSW-LMA 
scenario [l^l and background rate (B) from n C after the 
reduction with the technique outlined in the text. In Borexino 
the ratio between signal rate and raw background before any 
11 C rejection (B ) is 0.4; in KamLAND, S/B is 0.05; for a 
detector at SNOLab, S/Bo would be 36. 
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the two efficiencies C, V constrained to the value of e cor- 
responding to the chosen value of R. 

Further improvements of these figures might be 
achieved by the experiments by optimizing the set of cuts 
used to tag 11 C events with the double muon and neu- 
tron coincidence according to each detector's capabilities 
and performance. For example, the accurate reconstruc- 
tion of the track of the through-going muons might help 
reducing the dead mass-time by modifying the topology 
of the spatial cut proposed in this study, such as limit- 
ing the volume of the region excluded by the cuts to an 
intersection of a cylinder around the muon track and of 
a sphere centered around the neutron capture point, as 
proposed in [ie| . 

The technique offers as a by-product the possibility 
of determining the total number of 11 C decays in the 
detector. The residual background B can then be sta- 
tistically subtracted from the spectrum of the reduced 
data set. Let T be the total rate in the pep window after 
applying the muon+neutron tagging technique. Under 
the assumption that the internal and the external back- 
ground rates are negligible with respect to the residual 
background from 11 C and to the signal rates, S = T — B. 
A lower bound on the statistical error associated with the 
signal rate S extrapolated with the statistical subtraction 
of the energy spectra is obtained by propagating the er- 
rors in the formula just introduced: 



SS _ VN T + N B _ VN S + 2N B 
S ~ N s N s 

= J 1 + 2/R 0) 
VSM(l-D)to [> 

where Nt , Ns , and Nb are, respectively, the total number 
of events recorded, the total number of signal events and 
the total number of background events (all after applica- 
tion of the cuts to reduce the background from 11 C), to 
is the data taking time and M is the fiducial mass for pep 
and CNO neutrinos observation. The formula above can 
be used to obtain information concerning the statistical 
accuracy of the measurement of the signal rate. Borex- 
ino (KamLAND) with a 70-ton (300-ton) fiducial mass 
would achieve a statistical accuracy of 2.9% (3.4%) in 5 
years. 

VI. CONCLUDING REMARKS 

In this paper we presented a study of the production 
mechanism of 11 C nuclides in muon- induced showers. We 
identified the nuclear reactions relevant for the produc- 
tion of the nuclide in muon-induced showers in organic 
liquid scintillators. We performed an ab initio calculation 
of the production rates for each channel, then compared 
the calculated total production rate with available exper- 
imental data, obtaining a good agreement. We estimated 
that for 95% of the 11 C nuclides produced at least one 
neutron is emitted. 

A possible experiment located at SNOLab has a very 
low muon flux and hence a 11 C production rate which is 
negligible when compared with the expected rate from 
pep and CNO neutrinos. On the other hand, for both 
Borexino and KamLAND, cosmogenic 11 C is a signifi- 
cant background for the detection of pep and CNO solar 
neutrinos. 

We discussed a reduction technique for the 11 C events, 
based on one-by-one identification through the coinci- 
dence between a parent muon and the resulting neutron 
capture. We estimated that both Borexino and Kam- 
LAND could use the technique to increase the original 
signal (pep+CNO neutrinos) to background ( n C) ratio 
by a significant factor. Borexino can improve from a sig- 
nal/background ratio of 0.4 to one of 2 (4) while losing 
7% (58%) of the data to dead mass-time. KamLAND 
can improve its signal/background ratio from 0.05 to 0.2 
(0.3) while losing 12% (51%) of data to dead mass-time. 

The residual 11 C background can be statistically sub- 
tracted to determine the signal rate from neutrinos. We 
presented a formula providing a crude estimate of the 
statistical error associated with the signal rate deter- 
mined with this procedure. We estimated that both 
Borexino and KamLAND could measure the combined 
rate from pep and CNO neutrinos in the [0.8-1.3 MeV] 
window down to a 3% statistical accuracy in five years, 
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provided that internal and external background rates are 
kept within figures negligible with respect to the signal 
rate. 
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